home *** CD-ROM | disk | FTP | other *** search
/ Amiga Format CD 46 / Amiga Format CD46 (1999-10-20)(Future Publishing)(GB)[!][issue 1999-12].iso / -in_the_mag- / reader_requests / scilab / demos / lmitool / testnormopt < prev    next >
Text File  |  1999-09-16  |  457b  |  22 lines
  1. getf(SCI+'/demos/lmitool/normopt.sci','c')
  2. n=5;nu=2;ny=2;
  3. //rand('seed',0)
  4. A=rand(n,n);
  5. B=rand(n,nu);
  6. C=rand(ny,n);
  7. D=rand(nu,ny);
  8. //find H_inf norm of linear system A,B,C,D
  9. //   min gama 
  10. //   such that
  11. //   X-X'=0
  12. //  -[A*X+X*A',B,X*C';B',-gama*Ib,D';C*X,D,-gama*Ic]>0
  13. [X,gamaopt]=normopt(A,B,C,D)
  14. h_norm(syslin('c',A,B,C,D))
  15.  
  16.  
  17. getf(SCI+'/demos/lmitool/normoptd.sci','c')
  18. [X,g2opt]=normoptd(A,B,C,D);
  19. sqrt(g2opt)
  20. dhnorm(syslin('d',A,B,C,D),0.00001,1000)
  21.  
  22.